collapse_by_fe <- function(x, y, add_plus = TRUE) {
  if (add_plus) {
    return(paste0(paste(x, y, sep = " * "), " +\n    "))
  } else {
    return(paste(x, y, sep = " : "))
  }
}